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ABSTRACT 

Based on axi-symmetric hydrodynamical simulations and 3D reconstructions with Shape, we 
investigate the kinematic signatures of deviations from homologous ( "Hubble- type" ) outflows 
in some typical shapes of planetary nebulae. We find that, in most situations considered in 
our simulations, the deviations from a Hubble-type flow are significant and observable. The 
deviations are systematic and a simple parameterization of them considerably improves morpho- 
kinematical models of the simulations. We describe such extensions to a homologous expansion 
law that capture the global velocity structure of hydrodynamical axi-symmetric nebulae during 
their wind-blown phase. It is the size of the poloidal velocity component that strongly influences 
the shape of the position velocity diagrams that are obtained, not so much the variation of 
the radial component. The deviations increase with the degree of collimation of the nebula 
and they are stronger at intermediate latitudes. We describe potential deformations which these 
deviations might produce in 3D reconstructions that assume "Hubble-type" outflows. The general 
conclusion is that detailed morpho-kinematic observations and modeling of planetary nebulae can 
reveal whether a nebula is still in a hydrodynamically active stage (windy phase) or whether it 
has reached ballistic expansion. 

Subject headings: Planetary Nebulae, ISM: Outflows, Stars: Post-AGB, Stars: Mass loss, Hydrodynamics 



1. Introduction 

Most methods used for the reconstruction of 
three dimensional structures of planetary nebulae 
are based on imaging and kinematic data. Usually, 
they assume homologous expansion, i.e., purely ra- 
dial velocity with magnitude proportional to the 
distance from the central star (e.g. Sabbadin et 
al. 2006; Santander-Garcia et al. 2004; Hajian et 
al. 2007; Steffen & Lopez 2006; Morisset et al. 



2008; Dobrincic et al. 2008). The general justi- 
fication for adopting this so called "Hubble-type 
expansion" is that observed expansion speeds di- 
rectly correlate with distance from the center of 
the nebula (Wilson 1950, O'Connor et al. 2000, 
Meaburn et al. 2008). Furthermore, the struc- 
tures seen in position- velocity (P-V) diagrams are 
often very similar to the direct image (e.g. Mz-3: 
Santander-Garcia 2004; or NGC6302: Meaburn et 



1 



al. 2005). A technical reason for adopting this 
type of expansion law is that, it is the simplest di- 
rect mapping between the velocity of a gas parcel 
and its position along the line of sight. Still, this 
method requires an additional constraint to de- 
termine the constant of proportionality that links 
the Dopplcr- velocity to the position along the line 
of sight. This constraint may be an estimate of 
the extent of the object along the line of sight as 
deduced from partial symmetries or other argu- 
ments. 

The homologous expansion is a reasonable ap- 
proximation if the motion of the gas has become 
ballistic only a short time after the leaving the 
center of the nebula as compared to the current 
age of the nebula (for a discussion of this issue 
see e.g. Steffen & Lopez 2004). Observational 
and theoretical evidence is mounting that the hy- 
drodynamical events that determine the structure 
formation in many planetary nebulae are rather 
short-lived (e.g. Alcolea, Neri fc Bujarrabal 2007; 
Meaburn et al. 2005; Dennis et al. 2008; Akashi 
& Soker 2008). In order to decide whether the 
extended interacting wind scenario can be held up 
as it has been proposed almost two decades ago, it 
becomes important to look for kinematic evidence. 
During the action of the interacting winds, impor- 
tant deviations from a Hubble-type expansion law 
are to be expected in the velocity magnitude and 
direction. 

We study the kinematics of hydrodynamical 
simulations as a base to look for such deviations 
in observed nebulae. This study also aims to 
improve three-dimensional modeling that is be- 
ing carried out with photo-ionization codes (e.g. 
Morissct 2006; Ercolano, Barlow & Storey 2005) 
and morpho-kinematic reconstruction methods 
(e.g. Sabbadin et al. 2006; Stcffon & Lopez 2006). 

The morpho-kinematical modeling code Shape 
has been designed to allow the 3D modeling of very 
complex cmissivity structures and velocity fields 
(Steffen & Lopez 2006). Shape, and most current 
photo-ionization codes, do not actually calculate 
the structure and velocities from physical initial 
conditions. It is therefore important to have a 
physical guidance when deviations from a homol- 
ogous expansion are to be included in morpho- 
kinematical modeling. To provide improved veloc- 
ity fields we study the kinematics of some typical 
morphologies of planetary nebulae as created by 



two different mechanisms, purely hydrodynamics 
and magnetic confinement. 

As discussed by Li, Harrington & Borkowski 
(2002) detailed kinematic models are required to 
improve the accuracy of distance measurements to 
planetary nebulae from kinematic data. As these 
authors stress, deviations from homologous expan- 
sion should be taken into account to improve ac- 
curacy. Our paper aims at validating the Shape 
code for later application to the type of data pre- 
sented by Li, Harrington & Borkowski (2002) for 
BD-F30°3639. 

First, we apply purely hydrodynamical numeri- 
cal simulations to create elliptical and bipolar neb- 
ulae from an equatorial density enhancement in 
the environment. Highly elongated bipolar nebu- 
lae arc simulated with magncto-hydrodynamic cal- 
culations. The kinematics of the simulated data 
is analyzed and then reproduced as a morpho- 
kinematical model using Shape. We assess the dis- 
tortions that might occur in 3D reconstructions of 
such objects. In addition to providing template 
velocity fields, the direct comparison of the simu- 
lations with the reconstructed models using Shape 
allows us to validate the modeling capabilities of 
Shape itself. 

In Section 2 we describe the hydrodynamical 
model setup. Section 3 contains a description of 

the morpho-kinematical modeling technique and 
in Section 4 we discuss and summarize our results. 

2. Numerical simulations 

The four reference simulations have been per- 
formed using the magneto-hydrodynamic code 
ZEUS-3D (version 3.4), developed by M. L. Nor- 
man and the Laboratory for Computational Astro- 
physics. This is a finite-difi^erence, fully explicit, 
Eulerian code descended from the code described 
in Stone & Norman (1992). A method of charac- 
teristics is used to compute magnetic fields as de- 
scribed in Clarke (1996). The models include the 
Raymond & Smith (1977) cooling curve above 10'^ 
K. For temperatures below 10* K, the radiation 
cooling term is set to zero. Photoionization is 
not included for simplicity. The initial and mini- 
mum temperature allowed in all models is set to 
102 K. 

The computations are done in spherical coor- 
dinates (r, 9), with outflowing outer boundary 
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conditions and reflective conditions at the pole 
and equator. The models have grid resolutions of 
200 X 180 in r and 9, with a radial extent of 0.1 pc 
and 90°, respectively. The inner radial boundary 
is located at 2.5 x 10"'^ pc, i.e., we skip the inner 
2.5% of the computational mesh. 

The winds are set based on the rotating wind 
solutions given by Bjorkman & Cassinelli (1993), 
Ignace et al. (1996) and the approach to those 
equation given by Garcfa-Segura et al. (1999). 
These equations permit the introduction of the 
winds self-consistently. The magnetic field of an 
outflowing wind from a rotating star can be de- 
scribed by two components, and i?r (Chevalier 
& Luo 1994; Rozyczka & Franco 1996; Garci'a- 
Segura et al. 1999). The radial field component 
can be neglected since ^ while the toroidal 
component i?^ ~ r~^, and the field configuration 
obeys V • S = . 

Two useful dimensionless parameters contain- 
ing the basic information of the models are used. 
The first one, fi, is the velocity ratio of the stellar 
rotation to the critical rotation, and the second 
parameter, cr, is the ratio of the magnetic field en- 
ergy density to the kinetic energy density in the 
wind. As a initial condition, the AGB wind is 
set to fill up the whole computational mesh, while 
the fast wind is set only at the inner radial bound- 
ary, using equation (2) and (3) in Garcia-Segura 
et al.(1999) for both winds with their respective 
parameters. Note that equation (2) only sets the 
radial component of the wind velocity, since is 
negligible (at the inner boundary distance) due to 
angular momentum conservation. 

The four runs (Figure [T]) include, a bipolar (B) 
and an elliptical (E) case without magnetic fields, 
and a bipolar (Bm) and an elliptical (Em) case 
with magnetic fields on. The AGB slow winds 
have velocities of 10 km s~^ and M — 10^^ 
Mq yr~^ in all cases, while the fast winds have ve- 
locities of 1000 km s-i and M = 10"* Mq yr"! . 
Bipolar runs (B, Bm) have Vt = 0.98 in the AGB 
phases, while elliptical ones Vl ~ 0.6. Magnetic 
models have a = 0.01 only in the fast wind. AGB 
winds are set unmagnetized (see Garci'a-Segura et 
al. 1999 for further details). 
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Fig. 1. — Density distributions of model B at 2500 
yr (top-left), model E at 3500 yr (bottom- left) , 
model Bm at 1500 yr (top-right), and model Em 
at 2000 yr (bottom-right). 

3. Morpho-kinematical modeling with SHAPE 

In this study, the hydrodynamical simulations 
serve two purposes. First, we use them as observ- 
able objects, of which we know their structure, 
dynamics and kinematics. For a given time in the 
evolution, we analyze and visualize the projected 
images and position velocity (P-V) diagrams as if 
they where true observations. These images and 
P-V diagrams are generated with Shape after the 
data from the simulations have been imported in 
a suitable format (see below). 

The second purpose of using hydrodynami- 
cal simulations is to validate the functionality of 
Shape itself. For this purpose we reproduce the 
velocity field of the simulation and its structure in 
detail using the structure and velocity modeling 
tools of Shape . We then check whether the images 
and P-V diagrams are consistent with those of the 
simulation. 

We consider only the dense expanding shell that 
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has been swept up by the fast stellar wind. It 
is this shell, that is expected to be visible as the 
bright main shell in a planetary nebula. We as- 
sume that the nebula is fully ionized and that the 
emission measure, proportional to the square of 
the density, provides the relative brightness. 

In order to extract the region of the expanding 
shell we have selected cells from the hydrodynamic 
simulation according to their radial velocity and 
density contrast with neighboring cells. The se- 
lection is based on the fact that the shell is faster 
than the external medium (AGB-wind), but slower 
and denser than the active fast stellar wind. 

In Shape, each cell corresponds to a particle 
with properties of position, velocity and bright- 
ness. In order to convert the axi-symmetric two- 
dimensional hydro-data into a 3-D spatial particle 
structure, 18 particle copies for each extracted cell 
have been distributed randomly along the corre- 
sponding circle around the cylindrical symmetry 
axis. The velocity components of the particles 
have been adjusted according to the rotated po- 
sition. To visualize, analyze and further process 
the hydrodynamic simulations we import the sim- 
ulated particle data into Shape (Version 2.7). 

Since Version 1.0, Shape has an integrated tool 
for the creation of 3-D structures and their veloc- 
ity field. The simulated planetary nebula will be 
used to validate our modeling software, its method 
and accuracy. The hydrodynamic simulations that 
we are processing have been performed in spher- 
ical coordinates. Because of the cylindrical sym- 
metry, only the radial and the poloidal velocity 
components are non-cero. Shape has two different 
modeling modes for the velocity components. One 
is providing an analytic formula and the other is 
an editable curve, in which any number of linear 
segments can be used to describe the magnitude 
of each velocity component {vr,vg) as a function 
of the corresponding space coordinate (r, 9). 

The main shell of the nebula is modeled as a 
closed surface with cylindrical symmetry using the 
surface modeling tools in Shape . As an example. 
Figure [2] we show a visualization of the surface 
adjusted to run B of the simulations. The surface 
density of particles is constant within the devia- 
tions from the random particle statistics. About 
6 X 10"* particles have been used for each object. 
In order to model the total emission, the relative 
particle emissivity is given as a function of posi- 




Fig. 2. — The surface model mesh that has been 
adjusted to the extracted shell from run B. The 
dots that are plotted represent only 5 % of the to- 
tal particles that have been stochastically applied 
to sample the object. 

tion with a similar tool described for the velocity 
modeling (see above). 

For the rendering of the images, the emission 
is directly integrated along the line of sight. In 
this work we assume optically thin emission. The 
resulting raw image is then convolved with a gaus- 
sian point-spread-function that models the seeing 
conditions of an optical observation or the beam 
of a radiotelescope. 

Similarly, a P-V diagram is generated by dis- 
tributing the emission contained in a simulated 
spectrograph slit in a position-velocity image (ve- 
locity is horizontal and position along the slit is 
vertical). The raw image is then convolved with 
an ellipsoidal gaussian that has width and height 
corresponding to the velocity and spatial resolu- 
tion of the simulated instrument. 

4. Results and discussion 

In order to search for systematic tell-tale sig- 
natures which might reveal information about the 
physics of their initial formation process, we have 
investigated the kinematics of typical expanding 
planetary nebulae based on 4 hydrodynamic sim- 
ulations (Figured]). 
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Fig. 3. — The radial velocity component of the 
ellipsoidal and bipolar lobed hydrodynamical sim- 
ulations (E & B)is plotted against the distance 
from the center. The lines represent the piecewise 
linear representation of the velocity field in Shape. 
The dotted line is a simplified approximation with 
a single linear curve. The simulated data marked 
by (+) show those of model B and those marked 
as (x) are from model E. 
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Fig. 4. — The radial velocity component of the 
MHD-simulations is plotted against the distance 
from the center. For a description of the lines, see 
Figure ([3]). Here, two linear curve segments have 
been adopted for the simplified velocity field. The 
simulated data marked by (+) show those of model 
Bm and those marked as (x) arc from model Em. 

Figures [3] and 2] show the distribution of the 
radial velocities {vr) of the four simulations as a 
function of distance from the center. Figures [5] 
and [6] display the poloidal velocity (vg) compo- 
nent as a function of the angle from the polar 
axis of the simulations. For software validation 
the graphs have been approximated in Shape using 
piece-wise linear segments as shown by the con- 
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Fig. 5. — The poloidal velocity component of the 
ellipsoidal and bipolar lobed hydrodynamical sim- 
ulations is plotted against the angle in degrees 
from the symmetry axis (equator = 90°). For a 
description of the lines and symbols, see Figure 
([5]). Here, two linear curve segments have been 
adopted for the simplified velocity field. 



40 




-10 t . . . . . . . . 1 

10 20 30 40 50 60 70 80 90 
Angie (deg) 

Fig. 6. — The poloidal velocity component of the 
MHD coUimated simulations is plotted against the 
angle in degrees from the symmetry axis (equator 
= 90°). For a description of the lines and symbols, 
see Figure Here, three linear curve segments 
have been adopted for the simplified velocity field. 

tinuous lines with dots superimposed on the data 
from the simulations. The resulting P-V diagrams 
follow those of the simulation to within the scatter 
of the points from the simulation, validating the 
our Shape code. 

Figures ([7]) to (fTU)) show the rendered hydrody- 
namical simulations (columns "a"), their slit-less 
spectra or Hubble-like reconstructions (columns 
"b") and the reconstructed models with Shape 
(columns "c"). For the slit-less spectra, the simu- 
lated spectrograph slit was opened to include the 
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complete object. Note that a slit-less spectrum is 
equivalent to a 3D reconstruction that makes the 
assumption of a " Hubble- type" velocity distribu- 
tion, projected perpendicular to the line of sight. 
A small gaussian blur has been applied to simu- 
late observational seeing and reduce the stochastic 
particle-noise. We discuss the results for each sim- 
ulation separately. 

4.1. Non-MHD elliptical nebula 

Figures ([3]) and ^ show the statistics of the 
velocity component distribution in the hydrody- 
namical simulation of run E which produced an 
ellipsoidal nebula (Figure [7]). The radial velocity 
component increases approximately linearly with 
distance (Figure [3]). This might suggest that a 
homologous expansion is a good description of the 
velocity field. However, the plot of the poloidal 
velocity component versus position angle reveals 
that there are deviations of the velocity vector 
from the radial direction (Figured]). The devia- 
tions increase from the equator to poles until they 
reach approximately 7 km at 20 deg from the 
pole. Close to the poles the poloidal velocity de- 
creases again to zero. Where the poloidal compo- 
nent is largest, noticeable deformations will occur 
along the line of sight in a 3D reconstruction that 
assumes a homologous expansion law. 

The image, the slit-less spectra and the recon- 
structed models are compared in Figure ([7]). The 
shape of the slit-less spectrum at 90° (b) is more 
rectangular than the direct image (a). However, 
the reconstructed model (c) fits quite well. Al- 
though the deviations in (b) are small they are 
significant at intermediate to high latitudes. For 
the 0° and 90° viewing angles, the qualitative dif- 
ference is small and therefore might not influence 
significantly the physical interpretation of the ob- 
ject. However, at intermediate viewing angles the 
P-V diagram is more distorted and becomes point- 
symmetric. We consider it cautiously acceptable 
to directly reconstruct this type of objects from 
their spectral data assuming a homologous expan- 
sion law. 

Both, the radial and the poloidal velocity show 
small-scale fluctuations which are associated with 
local instabilities. In a reconstruction, unresolved 
fluctuations will result in an increased thickness of 
the reconstructed surface structure. 
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Fig. 7. — Column "a" shows the rendered images 
(emission measure) of the hydrodynamical simula- 
tion E at different inclination angles (0°,45°,90°) 
as seen from the observer's point of view. Column 
"b" contains the corresponding slit-less spectra of 
the simulation. It represents the view of the ob- 
ject as seen perpendicular to the line of sight in 
a reconstruction that is based on homologous ex- 
pansion. This image reveals the deformations ex- 
pected from such reconstructions. Despite of the 
cylindrical symmetry of the object, the shape of 
the " reconstruction" is quite different to the direct 
image of the simulation in column " a" . Column 
"c" shows the rendered image of our model as re- 
constructed with Shape (column " c" ) with correc- 
tions to the homologous expansion velocity field 
as seen from the same direction as column " b" . 

4.2. Non-MHD bipolar lobes 

The radial velocity distribution for the bipolar 
lobes from the purely hydrodynamic model B is 
shown in Figure ^ and the distribution of the 
poloidal velocity is given in Figure ([5]). As in the 
ellipsoidal case of run E, we find a nearly linear 
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Fig. 8. — The rendered images of the hydrody- 
namical simulation B (column "a"), slit-less spec- 
tra or homologous reconstruction (column "b") 
and the rendered image of the reconstructed model 
with Shape (column " c" ) . For further descriptions 
see Figure d?]). 

increase of the radial velocity component with dis- 
tance. Deviations from linearity are stronger far 
away and very close to the center. Also, the small- 
scale fluctuations increase at a certain distance. 
This distance marks the transition between two 
different dynamical regions in this bipolar nebula. 

The general trend of the poloidal velocity is 
similar to the one in the ellipsoidal simulation, a 
linear increase with angle up to a maximum value 
and then a decrease towards the pole. The van- 
ishing poloidal velocity at the pole and the equa- 
tor are enforced in these simulations due to the 
boundary conditions of reflection at the equator 
and polar axis. However, it is expected that this 
behavior is very closely sustained in real objects 
with approximate axi-symmetry. 

Column "b" in Figure ^ shows the slitless 
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Fig. 9. — The rendered images of the hydro- 
dynamical simulation Em (column "a"), slit-less 
spectra or homologous reconstruction (column 
"b") and the rendered image of the reconstructed 
model with Shape (column "c"). For further de- 
scriptions see Figure ([7]). 

spectra of simulation B for 3 different viewing an- 
gles. This corresponds to a 3D reconstruction that 
assumes homologous expansion as viewed along a 
direction perpendicular to the observer's line of 
sight. We note very significant differences to the 
observer's view in column "a". Especially, the 
"Hubble- type reconstructions" at 45° and 0° in- 
clination angles show strong deviations. As in the 
case of run E, the intermediate angles introduce a 
point-symmetry which is now even stronger. Col- 
umn "c" of Figure ([8|) shows the corresponding 
views of the axi-symmetric reconstructed model 
using Shape which applies the simple piece-wise 
linear velocity law (see below). 

In order to not only reproduce the images, but 
also the P-V diagrams with a model that improves 
the homologous expansion law with only a few ad- 
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Fig. 10. — The rendered images of the hydro- 
dynamical simulation Bm (column "a"), slit-less 
spectra or homologous reconstruction (column 
"b") and the rendered image of the reconstructed 
model with Shape (column "c"). For further de- 
scriptions see Figure ([7]). 

ditional parameters, we propose a simple piece- 
wise linear approximation to the velocity compo- 
nents as a function of position. 

For runs E and B, a two-parameter curve of 
two linear segments that start at zero at the pole 
and equator and meet at a point {9, ve) are suffi- 
cient. The point {9, vg) is a free parameter of the 
model. In addition, a linear relation of the radial 
velocity with distance could be used. The result of 
such a fit is shown in the P-V diagrams of Figure 
pn column E & B). The grey-scale P-V diagrams 
are from the hydrodynamical simulations with its 
original velocity information. The contours are 
from the same structural data, but with the new 
simple piece-wise linear velocity field. As an ex- 
ample, the corresponding velocity components of 
run B are plotted as dotted lines in Figures ([3|) & 
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Fig. 11. — Central slit P-V diagrams from the sim- 
ulations are represented in grey-scale. The super- 
imposed contours are from the same simulations, 
but with a simple piece-wise linear velocity dis- 
tributions applied (dashed lines in Figures ^ to 

m- 
m- 

We find that the agreement of this simple piece- 
wise linear velocity field with the simulation is very 
satisfactory. At any viewing angle deviations in 
the P-V diagram are only significant at the small 
scales of instabilities. Hence, for an elliptical or 
bipolar hydrodynamical nebula a very significant 
improvement in its morpho-kinematical modeling 
may be achieved by including only one and two lin- 
ear segments in the distribution of the radial and 
poloidal velocity components, respectively. This 
procedure retains the axi-symmetric structure of 
the object and, simultaneously, provides a good 
fit to the point-symmetric P-V diagrams. 

4.3. MHD collimated nebulae 

We now perform a similar analysis to the MHD- 
simulations. Figure shows the distribution of 
the radial velocity components of simulations Bm 
and Em as a function of distance from the cen- 
ter. The poloidal velocity component is plotted in 
Figure ([6|) . Both simulations show similar velocity 
fields, but with different positions and amplitudes 
of the dynamical regions. The radial velocity com- 
ponent has a nearly linear increase up to about 
half of the extent of the nebula. Then an abrupt 
change is noticeable with a small drop in the veloc- 
ity and further out a significant gradual increase 
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in the slope of the velocity graph. 

The poloidal velocity component is slowly grow- 
ing from the equator towards the pole and shows 
an accelerated increase from mid-latitudes on. It 
reaches a maximum near the tip of the outflow, 
after which it declines rapidly (as expected from 
the symmetry of the flow, see above). There are 
no strong small-scale variations in either velocity 
component as opposed to those found in the purely 
hydrodynamic simulations. 

Figures ([9]) and pO)l compare the projected im- 
ages rendered from the MHD-simulations (column 
"a") with the "3D Hubble-type reconstructions", 
i.e., the slit-less spectra (column "b"), and the 
reconstructed model using Shape (column "c"). 
The distortions observed in columns "b" are even 
stronger than in the non-MHD simulations. In- 
stead of being closed bipolar lobes, they appear to 
be open. In all cases of intermediate inclination 
angles, the axi-symmetry of the simulation is lost 
and replaced by point-symmetry. A physical in- 
terpretation of such artificial symmetry features, 
in 3D reconstructions that assume homologous ex- 
pansion, may lead to seriously erroneous results. 

In order to capture the main features of the 
velocity field in this type of objects in morpho- 
kinematical modeling, we propose to use at least 
two linear segments for the radial component and 
three segments for the poloidal velocity compo- 
nent. Such approximations for model Bm have 
been plotted in Figures (|4]) and ([6]) as dotted lines. 
See Figure pT|) to compare the P-V diagrams from 
the simplified Shape model with those of the sim- 
ulation. The deviations from the simulation are 
very small, hence, also for the MHD cases a con- 
siderable improvement of the velocity field can be 
obtained by the simple piece-wise linear velocity 
field proposed here. 

5. Proper motion patterns 

Additional kinematic information from an ex- 
panding object can be obtained from internal 
proper motion measurements. Such measurements 
have been obtained for a few nearby planetary 
nebulae (e.g. NGC6302: Meaburn et al. 2008; 
NGC7009: Rodriguez & Gomez 2007). Measure- 
ments of Doppler-shift and proper motion of a 
large number of features in a nebula (e.g. Li, Har- 
rington & Borkowski 2002) might reveal whether 




Fig. 12. — Proper motion vectors are shown for 
model B at 90 and 45 degrees viewing angle on 
the left and right, respectively. The top half has 
the velocity field as obtained from the simulation, 
while the bottom has a homologous expansion law. 

an object is hydrodynamically active or expanding 
homologously due to ballistic motion. In Figures 
(|12p and (|13p we compare the proper motion vec- 
tors of our simulations B and Bm with the cor- 
responding homologous motion patterns for two 
different viewing angles. The top halves of the im- 
ages represent the velocity pattern from the hydro- 
dynamical simulations, while in the bottom halves 
we replaced the original velocity field by a homol- 
ogous velocity field. Clearly, such measurements 
might reveal a lot of information about the dy- 
namical state of a planetary nebula. When the 
expansion is homologous, the proper motion pat- 
tern reflects this in purely radial vectors with mag- 
nitudes that increase linearly with distance. This 
pattern is basically independent of the structure 
and orientation of the object. 

However, in a dynamically active nebula or in 
a nebula with several kinematic sub-systems, the 
motion pattern will depend on the orientation and 
structure. Some of the proper motion vectors in 
Figures (|12p and (|13p almost reverse direction in 
the inclined objects (right side). Such a behavior 
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Fig. 13. — Proper motion vectors are shown for 
model Bm at 90 and 45 degrees viewing angle on 
the left and right, respectively. The top half has 
the velocity field as obtained from the simulation, 
while the bottom has a homologous expansion law. 
The top half includes part of the central coUimated 
outflow with very large velocity vectors. 

is found in BD+30°3639. As seen from the proper 
motion pattern in Li, Harrington & Borkowski 
(2002), BD+30°3639 is likely to still be a dynam- 
ically active nebula. The spectral and proper mo- 
tion data from Meaburn et al. (2008) indicate that 
NGC 6302 is expanding radially and, hence, pre- 
sumably ballistically. 

In general, because of the evolution of the fast 
stellar wind, it is expected that young planetary 
nebulae will be more likely to be active than 
older nebulae. It is therefore surprising that very 
young non-spherical nebulae, like, for instance, 
Minkowski's footprint (Ml-92, Alcolea, Neri & 
Bujarrabal 2007) show Doppler-shift patterns that 



are consistent with axi-symmetric homologous ex- 
pansion. It is important to systematically sur- 
vey the kinematic properties of planetary nebulae 
for evidence of deviations from homologous expan- 
sion. This will reveal information about the evo- 
lution of the fast stellar wind and the evolution 
of central stars in general. Hopefully, the analysis 
of large kinematical catalogs like those of Hajian 
et al. (2007) and Lopez et al. (2006) might be 
helpful in studying the transition phase between 
hydrodynamically active and ballistic expansion. 

6. Conclusions 

Based on axi-symmetric hydrodynamical simu- 
lations and 3D reconstructions with the software 
package Shape, we have investigated the kinematic 
signatures of deviations from "Hubble-type" out- 
flows in some typical shapes of planetary nebulae. 
The deviations are systematic and their descrip- 
tion as given in this paper may help to improve 
morpho-kineniatic and dynamic models of plane- 
tary nebulae. Excluding the small-scale variations 
from local instabilities, we find that the radial ve- 
locity as a function of distance can be well de- 
scribed by a linear increase for the cases of purely 
hydrodynamical ellipsoidal and bipolar- lobed neb- 
ulae. The poloidal component is, however, not 
negligable. Its parameterized approximation re- 
quires two linear segments, one increasing and the 
other decreasing with zero velocity at the pole 
and equator. We also find that the magneto- 
hydrodynamically generated cases require two lin- 
ear segments for the radial velocity component and 
three segments for the poloidal component. These 
simple extensions to a homologous expansion law 
capture the global velocity structure of hydrody- 
namically active axi-symmetric nebulae except for 
small-scale variations. 

We find that, in most situations considered in 
our simulations, the deviations from a homologous 
flow, both in magnitude and direction of the ve- 
locity vector, are significant and observable. The 
only case where only small deviations were found 
was the case of an ellipsoidal nebula. 

If significant deviations from a homologous ex- 
pansion model can be inferred, then it is very likely 
that it is still dynamically active in the sense that 
the fast wind is still significant or the inner region 
of the nebula is still over-pressurized with respect 
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to the exterior. In a future paper we will explore 
more precisely how the transition from the hydro- 
dynamically to the ballistic regimes occurs. 

It is the relative size of the poloidal velocity 
component that strongly influences the shape of 
the position velocity diagrams that are obtained, 
not so much the variation of the radial component 
with distance from the center. The deviations in- 
crease with the degree of collimation of the nebula. 
They are stronger at intermediate latitudes. There 
is a slow increase in deviations as we move from 
the equator to the poles, usually with a sudden 
increase at mid-latitudes. 

A good description of the velocity field in the 
MHD-simulations requires one more linear seg- 
ment than that of the non-MHD simulations. This 
would imply that this difference might provide ob- 
servational evidence for the presence or absence of 
dynamical active magnetic fields in the formation 
process of planetary nebulae. However, this aspect 
might require more detailed and specific studies. 

Summing up, detailed morpho-kinematic obser- 
vations and modeling of planetary nebulae can re- 
veal whether a nebula is still in a hydrodynami- 
cally active stage or whether it has reached ballis- 
tic expansion. 
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We acknowledge financial support from CONA- 
CYT grant 49447 and UNAM DGAPA-PAPIIT 
IN108506-2. Wc thank the anonymous referee for 
useful suggestions. 
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